Calorimetric Evidence of Multiband Superconductivity in Ba(Fe .925Coo.o75)2As 2 



F. Hardy/'BT. Wolf, 1 R. A. Fisher, 2 R. Eder, 1 P. Schweiss, 1 P. Adelmann, 1 H. v. Lohneysen, 1 ' 3 and C. Meingast 1 

1 Karlsruher Institut fur Technologie, Institut fur Festkorperphysik, 76021 Karlsruhe, Germany 
Lawrence Berkeley National Laboratory, Berkeley CA 94720, USA 
'Karlsruher Institut fur Technologie, Physikalisches Institut, 76128 Karlsruhe, Germany 

(Dated: October 26, 2009) 

We report on the determination of the electronic heat capacity of a slightly overdoped (x = 0.075) 
Ba(Fei_ a; Co a ;)2As2 single crystal with a T c of 21.4 K. Our analysis of the temperature dependence 
of the superconducting-state specific heat provides strong evidence for a two-band s-wave order 
parameter with gap amplitudes 2Ai (0)/ksT c = 1.9 and 2A2(0)/ksT c =4.4. Our result is consistent 
with the recently predicted s+_ order parameter [I. I. Mazin et al., Phys. Rev. Lett. 101, 057003 
(2008)]. 

PACS numbers: 74.25.Bt, 65.40.Ba, 74.20.Rp, 74.70.Dd 



The newly discovered iron arsenide family (FeAs) of- 
fers new possibilities for studying the interplay between 
superconductivity and magnetism^' 2 -' 3 ' As for many other 
materials, e.g., heavy fermions and cuprates, supercon- 
ductivity emerges in the vicinity of a magnetic insta- 
bility. The origin of the pairing interaction, as well as 
the gap symmetry remain unidentified in the pnictides. 
Theoretically, the particular topology of the Fermi sur- 
face with strong nesting features favor a multiband or- 
der parameter having either an s_| -wave or a rf-wave 

symmetryj^££ In cither case, a 7r-shift of the order- 
parameter phase is expected between different sheets 
of the Fermi surface. The identification of the gap 
symmetry is crucial, because it will shed light on the 
mechanism responsible for the condensation of Cooper 
pairs. Experimentally, solid evidence for a particular 
pairing state remains elusive, because several experimen- 
tal probes point to different conclusions. For instance, in 
the electron-doped 122 compound, Ba(Fei_ a: Co a ;)2As2, 
photoemission data 7 - (ARPES), and point-contact spec- 
troscopy 8 show two distinct nodeless gaps with large am- 
plitudes, while penetration-depth measurements- exhibit 
a power-law behavior reflecting the possible existence of 
nodes. Similar discrepancies are observed for hole-doped 
(Bai_ a; K a ;)Fe2As2 and the 1111 series i 10 ' 11 Some of these 
apparent contradictions may arise from the influence of 
the magnetic instability, which is expected to strongly al- 
ter the gap topology™ from impurity effects, or from ex- 
perimental difficulties like sample inhomogeneities or sur- 
face off-stoichiometry. Specific-heat measurements can 
provide an important measure of the bulk superconduct- 
ing properties; specifically, they can give valuable infor- 
mation about the possible existence of nodes in the en- 
ergy gap and, as previously shown for MgB2 , 13 ' 14 to the 
number of bands that contribute to the superconduct- 
ing condensate. Several specific-heat measurements have 
been reported for the Fe-pnictides, but the interpretation 
of the results has been impaired by substantial contri- 
butions from paramagnetic centers and/or an incorrect 
evaluation of the large phonon background ! 10 ' 15 ! 16 ' 17 

In this Letter, we present a detailed analysis of 
the electronic specific heat of a slightly overdoped 



Ba(Fei_ a; Co a ;)2As2 single crystal with x=0.075, i.e. at a 
doping level where the static spin-density wave (SDW) 
is no longer observed. The problem of the phonon 
background determination is overcome by measuring a 
strongly overdoped crystal with x=0.153, in which su- 
perconductivity is suppressed. Our analysis for the su- 
perconducting sample (x=0.075) gives strong evidence 
for two energy gaps, which implies that several sheets of 
the Fermi surface contribute to the formation of Cooper 
pairs. Additionally, we provide reliable values of the 
normal-state Sommerfeld coefficients j n =C e /T for sev- 
eral Co concentrations. 
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FIG. 1: Temperature dependence of the specific heat C of 
samples with x=0.075 and x=0.153. The dashed line repre- 
sents the lattice contribution, C; a i/T, derived from the spe- 
cific heat for x=0.153 (see text). The inset shows the low- 
temperature specific heat of both samples. 



Co-doped Bal22 single crystals were grown from 
sclf-flux in glassy carbon crucibles. Prereacted FeAs 
and CoAs powders were mixed with Ba, placed into the 
crucible, which then was sealed in an evacuated Si02 
ampoule. After heating to 650 °C and then to s» 1200 
°C with holding times of 5 hours, crystal growth took 



2 



place during cooling at a rate of « 1 °C/h. At 1000 
°C, the ampoule was tilted to decant the remaining 
liquid flux from the crystals and subsequently removed 
from the furnace. The composition of these samples was 
determined by energy dispersive x-ray spectroscopy to be 
x = 0.075 (± 0.005-0.01) and x = 0.153 (± 0.005-0.01), 
respectively. The specific heat was measured with the 
3 He option in a PPMS from Quantum Design. 

Figure [T] shows that the specific heat of both samples is 
dominated by the lattice contribution; the electronic part 
is only about rj 10 % of the total signal at T c . There- 
fore, it is impossible to obtain an accurate and unique 
description of the lattice background down to T = K 
using the harmonic-lattice approximation, i.e. by fitting 
the specific heat to an odd-power polynomial in a re- 
stricted range above T c . As will be demonstrated below, 
a much more reliable phonon specific heat, Gi a t (dashed 
line in FigfTJ) is obtained using the specific heat of the 
x=0.153 sample (for T>2 K) after subtraction of a con- 
stant electronic term j n = 15.8 mJ mol -1 K~ 2 . The elec- 
tronic term of this sample is obtained by fitting the data 
in the inset of Fig[T]with electronic and phononic terms, 
while ignoring the small peak at « 0.7 K, which may be 
due either to traces of remaining superconductivity or to 
the contribution from paramagnetic centers. It should be 
pointed out that the data of the superconducting sam- 
ple (x=0.075) contain a significant residual linear term 
of about 6 mJ mol -1 K~ 2 , which will be discussed later. 
No traces of the long range SDW were detected down 
to 0.5 K for either sample, in agreement with previous 
reportsi 17 ' 18 ! 19 
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FIG. 2: Temperature dependence of the electron specific 
heat, C e /T, of the superconducting sample (x=0.075). The 
dashed line represents the normal-state electron contribution, 
7„=23.8 mJ mol -1 K -2 . The dotted line is a residual normal- 
state-like contribution, 7 r =5.77 mJ mol -1 K -2 . 



Figure shows the temperature dependence of the 
electron specific heat C e (x=0.075) = C(x=0.075) - f s - 



Czat(x=0.153), with a scaling factor f s of 1.01. The fac- 
tor f s is introduced because it is not expected that the 
phonon specific heats of both samples are strictly iden- 
tical. The magnitude of f s was determined by enforc- 
ing entropy conservation, i.e. J^ c 7 n dT = J^ c C e /TdT. 
The small deviation of f s from unity demonstrates that 
the above procedure represents a very good method to 
determine the phonon background. Physically, this can 
be attributed to the fact that the substitution of Fe by 
Co does not substantially affect the lattice properties, as 
shown by recent inelastic x-ray scattering measurements 
and ab-initio calculations^ 



TABLE I: Critical temperature (T c ) and normal-state elec- 
tron specific heat (7 n ). Value for x=0 is taken from Ref.— 



X 


T C (K) 


7„(mJ mol" 1 K" 2 ) 





0.0 


5.3 


0.075 


21.4 


23.8 


0.153 


0.7 


15.8 



The superconducting transition at T c =21.4 K is re- 
markably sharp, indicating little inhomogeneity in the 
crystal. The normal-state electron contribution for 
x=0.075 is 7„ w 24 mJ mol -1 K~ 2 , which is in excellent 
agreement with LDA+DMFT calculations, which require 
7„ to be enhanced to 20-30 mJ mol -1 K~ 2 in order to 
explain mass renormalization by Hund coupling. 22 Our 
values for several Co concentrations, summarized in Ta- 
ble HI show that the disappearance of the SDW with Co- 
doping is accompanied by an increase of the electronic 
density of states (EDOS), compatible with a progressive 
closure of the SDW gap. In the overdoped region, on the 
other hand, j n and T c both decrease. Interestingly, j n of 
our superconducting sample is only about half as large as 
the value reported for K-doped 122 single crystals^ (« 
63 mJ mol" 1 K~ 2 ). Figure |1 illustrates that C e /T does 
not extrapolate to zero at T=0 but to a residual normal- 
state-like contribution 7 r =5.8 mJ mol -1 K -2 . Taken at 
face value, this would indicate that the sample has a su- 
perconducting fraction of « 75 %. Finite values of j r 
are a general feature of specific-heat measurements of 
electron-doped 122 iron arsenides i 10 i 23 For the cuprates, 
they have been attributed to an incomplete transition to 
the superconducting-state and volume fractions of nor- 
mal and superconducting material 7 r /7n and 1 - 7 r /7n, 
respectively. On this basis, the specific heat is the sum of 
separate contributions of the superconducting and nor- 
mal phases and consequently, the electronic specific heat 
can be normalized to one mole of superconducting mate- 
rial, C P 



^ es , defined by: 

C es = (C e 



IrT) 



In 



(1) 



However, recent specific-heat 23 and heat-transport 24 
measurements suggest that 7 r T is a consequence of pair 
breaking in electron-doped 122 pnictides, and not due to 
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the presence of normal material. Thus, if the s_| state 

is present^ 7 r T can be understood as arising from inter- 
band scattering, induced by in-plane disorder due to Co 
doping, which is pair-breaking for a sign-reversing order 
parameter^ In that case, the specific heat is, in prin- 
ciple, not the sum of contributions of broken pairs and 
the superconducting condensate. Nevertheless, in anal- 
ogy with the Na cobaltates^ C es (given by Eqfl]) can be 
expected to be a reasonable and useful approximation to 
the specific heat of one mole of superconducting material, 
and is used, in Fig. [31 for the purpose of comparison with 
several possible order parameters. Figure [3Ja) demon- 
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FIG. 3: a) The electron specific heat of the superconducting 
sample (x=0.075), normalized to 1 mol of superconducting 
condensate, compared with the specific heat of single-band s- 
wave (blue line) and d-wave (green line) order parameters, in 
the weak coupling limit, b) The electron specific heat of the 
superconducting sample (x=0.075), normalized to 1 mole of 
superconducting condensate. The red curve represents a two- 
gap fit. The blue and green curves are the partial specific-heat 
contributions of the two bands. 

strates that C es cannot be described by the specific heat 
of a single-band BCS s-wave superconductor, calculated 
in the weak-coupling limit (blue line). The agreement 
is very poor. As for MgB2r^ the positive curvature of 



C es for T/T c >0.6, where the BCS curve shows negative 
curvature, is indicative of strong coupling effects, but the 
observed discontinuity at T c , AC es /7„T c , which would 
be greater than the BCS value for a strong-coupled single- 
band superconductor, is close to weak-coupling value. In 
addition, C es is significantly larger than the BCS curve, 
for T/T c <0.4, again arguing against strong-coupling ef- 
fects. Figure[3Ja) also shows the specific heat of a single- 
band d-wave superconductor in the weak coupling limit 
(green line) ^ It is obvious that such a ^-dependent gap, 
even in a strong-coupling or a two-band scenario, can- 
not describe the observed low-temperature exponential 
behavior that can be inferred from the data. 

We therefore focus our discussion on the possibility of 
two energy gaps, using the phenomcnological two-band 
a-model, introduced by Bouquet et a L 14 ^ 27 It allows a 
fit of the specific heat from low temperatures up to T c 
and, as a result, gives reliable gap amplitudes that were 
shown to agree quantitatively with band calculations on 
MgB2 in particular ; 14 ' 28 and with Eliashberg equations 
in general^ In this fit, the specific heat is taken as the 
sum of contributions from two bands, which are calcu- 
lated independently assuming a BCS temperature de- 
pendence of the superconducting gaps. Two gap magni- 
tudes, at T=0 , are introduced as adjustable parameters, 
ai=Ai(0)/k B T c and a 2 =A 2 (0)/k s T c , together with a 
third quantity, ji/jn (*=1, 2), which measures the frac- 
tion of the total normal EDOS that the i-th band con- 
tributes to the superconducting condensate; 37 As shown 
in Fig[3jb), the two-band fit (with ai=0.95, «2=2.2 and 
7i/7n=0.33) accurately reproduces the specific heat over 
the entire temperature range. The smaller gap is about 
half the weak-coupling BCS value Abc , s(0) = 1.764, while 
the second gap is larger. These values are compara- 
ble with those derived from recent NMR (in the large 
scattering-rate limited) and /iSR-penetration-depth mea- 
surements^ (see Table UU, but differ appreciably, by at 
least a factor 1.5, from ARPES data! 7 - 



TABLE II: Gap ratios 2Ai(0)/k B T c , 2A 2 (0)/k s T c and 
weights 71/771 as determined by the two-gap model and by 
different techniques . 7 ' 30 ' 31 . 



Technique 


X 


2Ai(0)/k s T c 


2A 2 (0)/k s T c 


7i/7™ 


C (T) 


0.075 


1.9 


4.4 


0.33 


NMR 


0.070 


1.8 


7.2 


0.4 


fiSR 


0.070 


1.565 


3.768 


0.345 


ARPES 


0.075 


4.1 


6.4 





The temperature dependence of the superconducting- 
state specific heat, as well as the substantial resid- 
ual EDOS, is consistent with the predicted extended 
s-wave order parameter4££ However, the major gap 
develops around the Fermi-surface sheet that shows 
the largest EDOS while it is theoretically expected 
that A1/A2 oc yj N 2 /Ni, in the limit of pure inter- 
band pairing. 32 Thus, our results indicate that intra- 
band interactions are more important than expected. 
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Ba(Feo.925Coo.o75)2As2 also shares similar properties 
with NbSe2, another candidate for multiband super- 
conductivity, as illustrated by thermodynamic measure- 
ments^^ and ARPES spectral The specific heat 
of both NbSe2 and Ba(Feo.925Coo.o7s)2As2 show no 
sign of an incipient steep increase of C(T) below T c , 
which is the conspicuous signature of the small gap in 
MgB 2 , 13 This difference can be understood in terms 
of the gap anisotropy, A 2 /Ai, and the EDOS ra- 
tio, 72/71- This steep increase is particularly pro- 
nounced in MgB 2 because (i) A 2 /A! is about twice 
as large as in Ba(Feo.925Coo.o7s)2As2 and NbSe2, (ii) 
each gap gives an equal contribution, 72/71^1, to the 



specific heat of MgB2. In contrast, in both NbSe2 
and Ba(Feo.925Coo.o75)2As2 the major-gap contributions 
strongly dominate, with 72/71 roughly equal to 4 and 
2.3, respectively. 

In summary, a detailed analysis of the electronic spe- 
cific heat of Ba(Feo.925Coo.o75)2As2 provides strong ev- 
idence of a multigap order parameter, as observed for 
MgB2 and NbSe2- Our data are fit very well by a two- 
band s-wave model, and our results are, thus, in agree- 
ment with the predicted s_| pairing-state. Further, we 

derive a reliable phonon contribution that permits to ex- 
tract accurate values of the normal-state Sommcrfeld co- 
efficients. 
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